Density waves and supersolidity in rapidly rotating atomic Fermi gases 
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We study theoretically the low-temperature phases of a two-component atomic Fermi gas with 
attractive s-wave interactions under conditions of rapid rotation. We find that, in the extreme 
quantum limit, when all particles occupy the lowest Landau level, the normal state is unstable to 
the formation of "charge" density wave (CDW) order. At lower rotation rates, when many Landau 
levels are occupied, we show that the low-temperature phases can be supersolids, involving both 
CDW and superconducting order. 



The experimental achievement of condensation of pairs 
of atoms in two-component Fermi gases with resonant s- 
wave interactions [l|, [2|, y, |J] has allowed studies of inter- 
acting Fermi systems in regimes not accessible in solid- 
state systems: notably the transition region between 
weak and strong interactions (BCS to BEC crossover), 
and regimes of large density imbalance between the two 
species. The ability to rotate the gases, revealing a lat- 
tice of quantized vortices [f| , has provided an important 
diagnostic of superfluidity in these phase-coherent con- 
densates. 

A very interesting regime arises in atomic Fermi gases 
under conditions of rapid rotation (high vortex density). 
Noting the analogy between rotation and magnetic field 
in a superconductor, one might anticipate the BCS phase 
to revert to a normal state above a critical rotation fre- 
quency (analogous to H C 2 in superconductivity), as pre- 
dicted by BCS theory within a semiclassical approxima- 
tion [1, UM- Yet > 

going beyond this approximation to 
include Landau level (LL) structure, one finds that the 
normal phase can be unstable to ordered phases involving 
high-field superconductivity (SC) 0, "charge" density 
wave (CDW), or spin-density wave *(SDW) order 0. 

In this paper, we investigate the low-temperature 
phases of a two-component atomic Fermi gas with attrac- 
tive s-wave interactions under conditions of rapid rota- 
tion. The regime of interest for atomic gases differs sub- 
stantially from regimes studied in solid state systems: the 
rotation does not lead to any "Zeeman" splitting which 
might suppress high-field SC order; the short-range in- 
teractions allow density wave order to develop (this is 
suppressed in solid state systems by Coulomb interac- 
tions). We show that the low-temperature phases of an 
atomic Fermi gas with attractive interactions involve an 
interesting interplay between CDW and superconduct- 
ing phases. In the extreme quantum limit, when only 
the lowest Landau level (LLL) is occupied, we show that 
the system is unstable to CDW order along the rotation 
axis. At lower rotation rates, we show that CDW and 
SC can coexist, leading to "supersolid" behaviour. 

We study a rapidly rotating gas of two-species 
fcrmions, of equal densities, in the uniform limit: the 
number of vortices is assumed large, so the rotation fre- 
quency, Q, is close to the trap frequency, and the con- 



finement along the rotation axis is assumed weak. In the 
rotating frame, the Coriolis force mimics a magnetic field 
and leads to a Landau level structure with cyclotron fre- 
quency lo c — 2f2. The single particle states then have en- 
ergies e„ = (2n + l)hn + Ag^-, where v — (n, x, k) stands 
for the LL-index n, the momentum in the Landau gauge 
x [24|, and the wavevector along the rotation axis k. For 
a non- interacting gas with Fermi energy €f , the n th Lan- 
dau level has a ID Fermi surface with Fermi momentum 
fikpn = [2m(eF — (2n + l)fi.f2)] 1//2 and kinetic energy rel- 
ative to the bottom of the band ep n = % 2 k 2 7n /2m. We 
describe the instabilities of these Fermi surfaces arising 
from weak interactions. (We focus on results for attrac- 
tive interactions, but also report on the repulsive case.) 

First, we analyze the effect of rotation on the SC phase, 
applying BCS theory in the presence of Landau level 
structure [H, 0, H, EH H, 13 1. For contact interactions, 
the gap equation requires regularisation at high energies. 
For solid state systems, the Debye frequency provides 
a natural cut-off for phonon-mediated attractive interac- 
tions. In a cold atomic gas, a natural regularization arises 
from the (small) lengthscale of the interparticle forces. 
Using a two-channel model for the Feshbach interaction, 
this lengthscale enters as the size of the "closed channel" 
boson (see e.g. [1, EH, EH) and can be taken to zero with 
the introduction of appropriate counterterms. Following 
Ref. [3], the parameters of the model are the boson en- 
ergy cb = 2hfl + S+C and the coupling aS vv i between a 
closed channel boson and fermions with quantum num- 
bers v and v' . Here, S is the physical detuning of the 
bosons and C a counterterm which is set to cancel the 
boson self energy £(w -> 0) = a 2 \S VV > \ 2 /(e„ +e v >), 
such that the model reproduces the scattering properties 
at low energy and f2 — > [3][2H|- The physical scatter- 
ing parameters are related via —a 2 /8 = Airh 2 a s / m = g. 
Treating the ensuing two-channel Hamiltonian within 
mean field, and assuming a wide Feshbach resonance, 
yields the linearized gap equation [idj ] 
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FIG. 1: (color online) The critical temperature calculated 
within BCS mean-field, for rotation frequency £1 and chemical 
potential /i. Arrows indicate critical frequencies Q C 2 to the 
left of which T c vanishes in the semiclassical approximation 
0|- Dashed lines show T c of the CDW state in the LLL as 
obtained in the parquet approximation. 



length £ = (h/2mfl) 1 / 2 . The solutions to |T]) determine 
the critical temperature T c for superconductivity. 

Within a semiclassical approximation to Qj_ T c van- 
ishes for Ml > A 2 /M (A is the zero field gap)[S 0, Q [1| • 
The full gap equation (1) admits solutions even in this 
regime. Then, when T c is small, the dominant contribu- 
tions arise from integrating the 'diagonal' terms (n = n') 
[l2T ] , which diverge logarithmically at low T for occupied 
LL's. Provided k B T c < [fj, - TiCl(2n max + 1)], the off- 
diagonal terms (n ^ n') can be neglected, and one finds 
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The critical temperature ([2|) is a strongly oscillating func- 
tion of [i/h£l, with a peak each time a LL depopulates 
and G(rj) diverges. The sharp peaks predicted by (|2|) 
are rounded in a full of solution of ([1]) which is required 
for strong-coupling. The evolution from weak to strong 
coupling is shown in Fig. [TJ which we have computed by 
solving (p} using a numerical root-finding routine. 

Consistent with previous studies of BCS theory in 
solid state systems @, 0, S El 12, 13 1 we find that 



LL quantization leads to a stable SC state at any value 
of the field 9]. For \a s \kp ^ 1 the critical tempera- 
ture has a minimum value T™ m (a s fci?). For temperatures 
T > T™ m (a s fci?), mean-field theory results predict a se- 
ries of reentrant SC to normal transitions as the rotation 
rate increases. Our results differ from those presented 
in Ref. the (reentrant) superconductivity at rapid 

rotation was not found in that work; the critical rotation 
frequency has an important temperature-dependence, es- 
pecially for strong coupling. 

The superconducting phase competes with other or- 
dered phases. To determine the nature of the groundstate 
one must work beyond mean-field theory. We analyze the 
competition between SC and other ordered phases within 
a single LL, for example when all particles occupy the 
LLL. Owing to the quasi-lD dispersion within this LL, 
any response function that connects opposite sides of the 
Fermi surface diverges at low temperatures as a power of 
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Identical divergences occur in both particle-particle (p- 
p) and particle- hole (p-h) diagrams [1 €31 ] , and in diagrams 
of higher orders. The resulting ensemble of "parquet" 
diagrams, obtained by mutual insertion of p-p and p-h- 
blocks into one another [171 ]. is most easily analyzed in 
terms of a renormalization group (RG) approach [lij ]. 
This scheme has been applied to spinless electrons in a 
magnetic field [lij ]. We generalise this approach to a 
iwo-component rotating atomic Fermi gas, and consider 
particles in the n th LL. Particles at the two Fermi points 
k = ±fci?n are represented by separate fermionic field 
operators a^) and b">. States at the first Fermi point 
are expanded in terms of the LL wavefunctions 24 1 



Vw(*, Y, Z) = Af n H n (X - x y*r+(x~ x f/2 +l kz^ (5) 

with lengths measured in units of £q, the Hermite poly- 
nomials H n and normalization J\f n = (L v L z tt^ 2 n nl)~ 2 . 
At the other Fermi point we use the transformed basis 

i> nyk (X, Y,Z) = ^-J2 (X, Y,Z). (6) 



For weak coupling, the kinetic energy can be linearized 
around the Fermi points, ±kp n . The (logarithmically 
divergent) part of the contact interaction (amplitude g) 
describing scattering between opposite Fermi surfaces is 
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The dependence on the LL index n arises only in 
the form of the bare interaction vertex, 7q™' ) (r) = 
e~5 r2 [i„(r 2 /2)] 2 , where we introduce r = (x, y) and L n 
are the Laguerre polynomials. The interaction Hi can be 
viewed as two distinct vertices according to the way spin 
is conserved, and denoted 7^2 in the usual notations for 
quasi- ID systems From ([7]), these vertices have the 
initial conditions 



7i™2( r )l?=o = sgn( 5 )7^ n) (r). 
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Renormalisation of the vertices 71,2 leads to corrections 
that can be expressed as a power in £ Q 17]. The one- 
loop RG equations can be obtained by adapting the ap- 
proach of Ref. [20] to include the LL structure. We find 



—— = -2 71 * 71 + 2 71 * 72 
d( 
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-2 7i®72 (9) 
-71 ® 71 - 72 <8> 72 (10) 



where the operations * and ® arise in p-h and p-p loops, 
respectively, and are defined by 
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The phase factor in (fT2"]) is a consequence of the LL struc- 
ture. 

We have solved the RG equations (|9I10|) with initial 
conditions © for arbitrary Landau level index n, using 
a standard numerical routine with 7i j 2(l r |) discretized 
uniformly in |r|. (The initial conditions for 7^2 are ra- 
dially symmetric, and this symmetry is preserved by the 
RG equations.) To identify instabilities, we calculate the 
renormalization of the response functions [23]. The RG 
equations for the triangular vertices T in the (singlet) 
SC, charge- and spin-density wave (SDW) channels are 
given in our case by 



d^Tssc = (-71 - 72) <S> Tssc 
<%^cdw = (-271 + 72) * Tcnw 

= 72 * ^SDW- 
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Initial conditions for the triangular vertices can be cho- 
sen as %\^ =0 = <5(r), such that all Fourier components 
are non-zero. We find the smallest value, £ c , at which a 
susceptibility diverges: this indicates a transition into an 
ordered phase at a critical temperature [see (j3J)] 
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In contrast to the full RG equations, the simplified 
equations describing only p-h ladders can be solved an- 
alytically, and provide a useful reference point for our 
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TABLE I: With dynamics restricted to a single Landau-level, 
the analysis of the parquet diagrams reveals a CDW instabil- 
ity for attractive interactions g < 0, and a SDW instability 
for g > 0. While the CDW is enhanced by scattering in the 
p-p-channel, SDW order is weakened. As n — > 00 the critical 
values, £ c , converge to the result for p-h-ladders, £ c = (2tt)~ 1 . 



numerical evaluation. The solution for the p-h ladder 
discussed in Ref. [HI] can be generalized to arbitrary LL 
index n and yields a transition at a critical tempera- 
ture which is independent of n and the sign of g, with 
£ c = (27r) _1 . [For g < (g > 0) the transition is to 
a CDW (SDW).] For p-p ladders, the problem can be 
solved analytically for n = 0, where the SC instability 
occurs for attractive interactions also at £ c = (27r) _1 . By 
restricting the SC gap equation in the presence of a mag- 
netic field to a single LL (see above and [13]), one can 
infer f c (n) = (2 n ra!) 2 /[27r(2n)!], showing that SC order 
becomes weak as n — > 00. These analytic results are re- 
produced by our numerical approach, when restricted to 
include p-p or p-h diagrams only. Note that for n = the 
CDW and SC instabilities have the same critical temper- 
ature. Thus, mean-field theory cannot determine which 
of these states will form the low-temperature phase. 

Our solution of the full RG equations (|9ll0p shows 
that, for attractive interactions, CDW order is the dom- 
inant instability for all LLs. The critical temperature 
(|16p depends on the LL index, with exponents summa- 
rized in Table [U Thus, for the LLL n = 0, the competi- 
tion between the identical instabilities in p-p and p-h 
channels [both at £ c = (27r) -1 ] is decided to the ad- 
vantage of CDW order. The order parameter diverges 
most strongly at zero in-plane momentum, so the density- 
waves are aligned with the rotation axis. Thus, the 
CDW phase in the n th LL involves a modulation of 
the particle density along the rotation axis, with pe- 
riod A^ DW = n/kp n . Within one period of the den- 
sity wave the effective 2D particle density (in that LL) 
is n.2d,n = l/( 7r ^o)' snch that this LL is fully occupied 
(its filling factor is v n = n 2 d.ri27r^Q = 2). Thus the CDW 
phase is fully gapped. In the extreme quantum limit, 
when HQ < ep < 3hQ, fc^o — 7r 2 n£§, where n is the 3D 
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particle density, so the period is Aq DW = l/( 7r "^o)- I 11 
Fig. Q] we show the transition temperature into this CDW 
in the LLL (dashed lines). (For repulsive interactions, we 
find that SDW order is dominant for all n. See Table |U) 

Our results show that CDW order always prevails for 
attractive contact interactions when dynamics are re- 
stricted to a single LL. However, at low rotation rates, 
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the groundstate is the BCS superconducting state (with 
dilute vortices). How does one reconcile these conclu- 
sions? The answer lies in the coupling between LLs. 
Since the periods of the CDWs, A^ DW , differ between 
LLs, we find that the CDW does not gain from inter-LL 
couplings: there are CDW instabilities at the tempera- 
tures set by our calculations for individual LLs, Table fl] 
On the other hand, a SC state can benefit from coher- 
ence between LL's, as the Cooper pairs all have the same 
(zero) momentum. Thus, although SC within a single LL 
is less relevant than CDW, the "Josephson" coupling be- 
tween LLs can stabilise a collective SC state. That said, 
as the topmost LL, n max , depopulates our results show 
that the CDW instability in this LL can occur at a higher 
temperature than the SC state of the entire system. In 
this case, the first instability (as T is reduced) is to a 



CDW in the Landau level 



and one expects a sec- 



ond instability, at lower T, to a SC state formed from the 
other Landau levels. (The loss of the highest LL from the 
SC makes little difference to its condensation energy.) In 
this way, we predict a supersolid groundstate, involving 
both CDW of the topmost LL and SC order in the lower 
LLs. Ultimately, at sufficiently high rotation rate (or low 
particle density), when all particles occupy the LLL, the 
groundstate is a CDW without superconducting order. 

A striking consequence of our results is that for a 
rapidly rotating atomic Fermi gas, there should appear 
spontaneous density wave order, with a period A^ DW that 
grows as the particle density in the topmost Landau level 
decreases. This can be a long lengthscale, so could be 
measured in experiment directly by in situ absorption. 
Clearly, the observation of the density waves requires a 



CDW 



For t\\ < A 



CDW 



trap with oscillator length £^ > A 
there will be a single period of the wave, leading to a 
quasi-2D regime with 2D particle density in this LL equal 
to U2d,n = l/( 7r ^o)- This (incompressible) filled LL will 
appear as a step in the transverse density profile, as mea- 
sured in-situ or in an expansion measurement 
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The results that we have presented are accurate far 
from the resonance on the BCS side, where interactions 
are weak. We expect the qualitative behaviour to survive 
as the resonance is approached. While the detailed en- 
ergetics of both phases cannot be relied upon for strong 
coupling, we find that SC is stabilized relative to CDW 
order for chemical potentials above the LLL as the cou- 
pling increases. Presumably, this leads to the suppres- 
sion of CDW states in any but the lowest LL as one 
approaches the resonance. Furthermore, we note that 
the density-wave state(s) we find on the BCS side of the 
resonance cannot evolve smoothly to the BEC side. A 
CDW of atoms, with u aiom = 2 per period, could evolve, 
to retain the same period, into a CDW of tightly bound 
molecules with v mo \ = 1/2 per period [22| . However, 
there must be a phase transition separating these two 
states, owing to the different edge structures of the phases 
Thus, in contrast to the SC phase at low rota- 
tion rate, in the extreme quantum limit (at high rotation 



rate) tuning the interactions across the Feshbach reso- 
nance must involve a phase transition. 
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withx = (x-x')/2, N = n+n', and / dX / (2-k)\c x \ 2 = 1. 
[26] At T = 0, the kinetic energy of the superfluid flow in the 
vortex lattice (in the rotating frame) is ~ HQ, per particle, 
while the condensation energy is ~ A 2 //i per particle. 



